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The recent realization of a "Levy glass" (a three-dimensional optical material with a Levy dis- 
tribution of scattering lengths) has motivated us to analyze its one-dimensional analogue: A linear 
chain of barriers with independent spacings s that are Levy distributed: p(s) oc s -1-< * for s — » oo. 
The average spacing diverges for < a < 1. A random walk along such a sparse chain is not a 
Levy walk because of the strong correlations of subsequent step sizes. We calculate all moments 
of conductance (or transmission), in the regime of incoherent sequential tunneling through the bar- 
riers. The average transmission from one barrier to a point at a distance L scales as L~ a InL for 
< a < 1. The corresponding electronic shot noise has a Fano factor (oc average noise power / 
average conductance) that approaches 1/3 very slowly, with 1/lnL corrections. 

PACS numbers: 42.25.Dd, 05.40.Fb, 42.68.Ay, 73.50.Td 



I. INTRODUCTION 

In a recent publication Barthelemy, Bertolotti, and 
Wiersma have reported on the fabrication of an unusual 
random optical medium which they have called a Levy 
glass. It consists of a random packing of glass micro- 
spheres having a Levy distribution of diameters. The 
space between the spheres is filled with strongly scat- 
tering nanoparticles. A photon trajectory therefore con- 
sists of ballistic segments of length s through spherical re- 
gions, connected by isotropic scattering events. A Levy 
distribution is characterized by a slowly decaying tail, 
p(s) oc l/s 1+Q for s — ► oo, with < a < 2, such that the 
second moment (and for a < 1 also the first moment) di- 
verges. The transmission of light through the Levyglass 
was analyzed [l| in terms of a Levy walk 0, d, 0] for 
photons. 

Because the randomness in the Levy glass is frozen 
in time ("quenched" disorder), correlations exist be- 
tween subsequent scattering events. Backscattering af- 
ter a large step is likely to result in another large step. 
This is different from a Levy walk, where subsequent 
steps are independently drawn from the Levy distribu- 
tion ("annealed" disorder). Numerical [5| and analytical 
Q theories indicate that the difference between quenched 
and annealed disorder can be captured (at least approx- 
imately) by a renormalization of the Levy walk expo- 
nent — from the annealed value a to the quenched value 
a' = a + (2/rf) max(0, a — d) in d dimensions. Qual- 
itatively speaking, the correlations in a Levy glass slow 
down the diffusion relative to what is expected for a Levy 
walk, and the effect is the stronger the lower the dimen- 
sion. 

To analyze the effect of such correlations in a quan- 
titative manner, we consider in this paper the one- 
dimensional analogue of a Levy glass, which is a lin- 
ear chain of barriers with independently Levy distributed 
spacings s. Such a system might be produced artificially, 
along the lines of Ref. [l|, or it might arise naturally 
in a porous medium 0] or in a nanowire 0]. Earlier 



studies of this system 0, [Io|, HH, EH have compared 
the dynamical properties with those of a Levy walk. In 
particular, Barkai, Fleurov, and Klafter [llj found a su- 
perdiffusive mean-square displacement as a function of 
time [(x 2 (t)) oc V with 7 > 1] — reminiscent of a Levy 
walk (where 7 = 3 — a). No precise correspondence to 
a Levy walk is to be expected in one dimension, because 
subsequent step lengths are highly correlated: Backscat- 
tering after a step of length s to the right results in the 
same step length s to the left. 

The simplicity of one-dimensional dynamics allows for 
an exact solution of the static transmission statistics, 
without having to assume a Levy walk. We present 
such a calculation here, and find significant differences 
with the L~ a l 2 scaling of the average transmission ex- 
pected [13, [H, [lj| for a Levy walk (annealed disor- 
der) through a system of length L. If the length of 
the system is measured from the first barrier, we find 
for the case of quenched disorder an average transmis- 
sion (T) oc L~ a bxL for < a < 1 and (T) oc L~ x for 
a > 1. Note that the nonalgebraic length dependence 
for < a < 1 goes beyond what can be captured by a 
renormalization of a. 

In the electronic context the average conductance (G) 
is proportional to (T), in view of the Landauer formula. 
In that context it is also of interest to study the shot 
noise power 5, which quantifies the time dependent fluc- 
tuations of the current due to the granularity of the elec- 
tron charge. We calculate the Fano factor F oc (S)/(G), 
and find that F approaches the value 1/3 characteristic 
of normal diffusion fl7l . Il8j with increasing L — but with 
relatively large corrections that decay only as 1/ In L for 
< a < 1. 



II. FORMULATION OF THE PROBLEM 

We consider a linear chain of tunnel barriers, see Fig. [I] 
with a distribution of spacings p(s) that decays for large 
s as 1/ s a+l . A normalizable distribution requires a > 0. 
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FIG. 1: Linear chain of randomly spaced tunnel barriers. We 
study the statistics of conductance (or transmission) over a 
length L for a Levy distribution of spacings p(s). 



For < a < 1 the mean spacing is infinite. We take 
for each barrier the same mode-independent transmission 
probability T <C 1 (no ballistic transmission) . The corre- 
sponding tunnel resistance is r = (h/e 2 )(NT)~ 1 , with iV 
the number of transverse modes. In the electronic con- 
text we require r <C h/e 2 , so that the Coulomb blockade 
of single-electron tunneling can be ignored. 

We work in the regime of incoherent sequential tunnel- 
ing (no resonant tunneling). This regime can be reached 
for N ^> 1 as a result of intcrmode scattering, or it can be 
reached even for small iV as a result of a short phase co- 
herence length. For sequential tunneling the resistance R 
of n barriers in series is just the series resistance nr [cor- 
responding to a transmission probability T = (nr) -1 ]. 
We measure this resistance 



R(L) = rJ26{x n )0(L-x n ) 



(2.1) 



between one contact at x — and a second contact at 
x = L > 0. The numbers x n indicate the coordinates of 
the tunnel barriers and 6{x) is the step function [6(x) = 1 
if x > and 0(x) = if x < 0]. 

Without further restrictions the statistics of the con- 
ductance would be dominated by ballistic realizations, 
that have not a single tunnel barrier in the interval (0, L). 
The reason, discussed in Ref. [ll| , is that the average dis- 
tance between a randomly chosen point along the chain 
and the nearest tunnel barrier diverges for any < a < 2 
(so even if the mean spacing between the barriers is fi- 
nite). To eliminate ballistic transmission, we assume that 
one tunnel barrier is kept fixed at xq = + . (This bar- 
rier thus contributes r to the resistance.) If we order the 
coordinates such that x n < x n+ i, we have 



We seek the scaling with L in the limit L — > oo of the 
negative moments (R(L) P ) (p = —1, —2, —3, . . .) of the 
resistance. This information will give us the scaling of 
the positive moments of the conductance G = R^ 1 and 
transmission T — (h/Ne 2 )R . It will also give us the 
average of the shot noise power S, which for an arbitrary 
number of identical tunnel barriers in series is determined 
by the formula [19( 



S=-e\V\r- 1 l(R/r)- 1 +2(R/r)-% (2.3) 



where V is the applied voltage. From (S) and (G) we 
obtain the Fano factor F, defined by 



F = 



(S) 



2e\V\(G) 



(2.4) 



III. ARBITRARY MOMENTS 



The general expression for moments of the resistance 



is 



R(L) 



J2e(x n )e(L- Xn ). 



(2.2) 



(B(^> = H , ni + £8(i n )9(I-i B )J V (3.1) 

where the brackets (•••} indicate the average over the 
spacings, 

°° />oo 

(•'•) = H/ dx n p(x n -x n -\) ■■■ , (3.2) 

n=1 J-oo 

with the definitions xq — and p(s) — for s < 0. We 
work out the average, 

(RiLf) = rPj2 nP [Ii / dsiPisi)) 

x 9 [J2 Si - l\0 I L - ^ sA . (3.3) 



It is more convenient to evaluate the derivative with 
respect to L of Eq. (|3.3p . which takes the form of a mul- 
tiple convolution of the spacing distribution [2fj| . 
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In terms of the Fourier (or Laplace) transform 
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the series (|3.4[) can be summed up, 
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The function Li(cc) is the polylogarithm. The imaginary 
infinitesimal i0 + added to £ regularizes the singularity of 
the integrand at £ = 0. For negative p this singularity is 
integrable, and the integral (|3.6p may be rewritten as an 
integral over the positive real axis, 



(3.7) 



IV. SCALING WITH LENGTH 
A. Asymptotic expansions 

In the limit L — > oo the integral over £ in Eq. (|3.7[) is 
governed by the £ — » limit of the Fourier transformed 
spacing distribution. Because p(s) is normalized to unity 
one has /(0) = 1, while the large-s scaling p(s) cx l/s a+1 
implies 

lim f(0 = ( ; + C " (S °^ )Q ; ° < a < X ' (4.1) 
^o J ^' \ 1 + «s£ + c Q (s £) Q , 1 < a < 2. v 7 

The characteristic length sq > 0, the mean spacing s, as 
well as the numerical coefficient c a are determined by the 
specific form of the spacing distribution. 

The limiting behavior of the polylogarithm is governed 

by 

Lii(l + e) = -In(-e), (4.2) 
limLi 2 (l + e) =C(2)-eln(-e), (4.3) 

limLi n (l + e) = C(n)+C(n-l)£, n = 3,4,... (4.4) 

In combination with Eq. (|4.1[) we find, for < a < 1, the 
following expansions of the integrand in Eq. (|3.7[) : 
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limi--^Li_ p (/)= -C(-p)c Q ( So O Q , 
p=-2,-3... 



(4.5) 
(4.6) 



For 1 < a < 2 we should replace c Q (so£) Q by is£ 

Ca(soO Q - 



B. Results 



We substitute the expansions (|4. 5|) and (|4.6|) into Eq. 
(|3.7p . and obtain the large- L scaling of the moments of 
conductance with the help of the following Fourier inte- 



TABLE I: Scaling with L of moments of conductance (or, 
equivalently, transmission) . 
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grals (L > 0, a > -1): 
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Here je is Euler's constant and H a is the harmonic num- 
ber. The resulting scaling laws are listed in Table Q] 
Two physical consequences of these scaling laws are: 

• The Fano factor (|2.4|) approaches 1/3 in the limit 
L — > oo, regardless of the value of a, but for < 
a < 1 the approach is very slow: F — 1/3 oc 1/ In L. 
For 1 < a < 2 the approach is faster but still sub- 
linear, F — 1/3 oc l/i"" 1 . 

• The root-mean-squarc fluctuations rms G = 
a/ (G 2 ) — (G) 2 of the conductance become much 
larger than the average conductance for large L, 
scaling as rms G/{G) oc L a / 2 /\nL for < a < 1 
and as rms G/(G) oc L 1 -"/ 2 for 1< a < 2. 



V. NUMERICAL TEST 

To test the scaling derived in the previous sections, 
in particular to see how rapidly the asymptotic In- 
dependence is reached with increasing L, we have nu- 
merically generated a large number of random chains of 
tunnel barriers and calculated moments of conductance 
and the Fano factor from Eqs. (|2~2"]) - (|2~4"1, . 

For the spacing distribution in this numerical calcula- 
tion we took the Levy stable distribution 2l| for a = 1/2, 



Pi/aOO = ( S o/2 7 r) 1 /2 s -3/2 e - S o/2^ 



(5.1) 



Its Fourier transform is 



A/a(0 = exp(-V-2is => c 1/2 = i - 1. (5.2) 
Inserting the numerical coefficients, the large-L scaling 



4 




0.01 



10 100 1000 

L/s 



10000 



FIG. 2: Scaling of the average conductance (bottom panel), 
the variance of the conductance (middle panel) , and the Fano 
factor (top panel), for a chain of tunnel barriers with spacings 
distributed according to the a = 1/2 Levy stable distribution 
(|5.1|) . The data points are calculated numerically, by averag- 
ing over a large number of random chains of tunnel barriers. 
The solid curves are the analytical results (|5.3[) - (|5.6[) of the 
asymptotic analysis in the L — > oo limit. 



of conductance moments for the distribution (|5.1[) is 
1 



lim (G) = -(2^L/ So r 1/2 [ln(2L/ So ) +1e] 



1 



lim {(?) = 2C(p)-(27tL/ So )- 1/2 , p > 2. 
L^oo r p 



(5.3) 
(5.4) 



The resulting scaling of the conductance fluctuations and 
Fano factor is 



F 



rms G 

1 

3 4 



_{G^ 
(4/3)C(3) 



- 1 



(^/3)(2^L/ So ) 1 / 2 
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ln(2L/s ) + IE ' 



1, 

(5.5) 
(5.6) 



In Fig. [2] we compare these analytical large- L formu- 
las with the numerical data. The average conductance 



converges quite rapidly to the scaling (|5.3|) . while the 
convergence for higher moments (which determine the 
conductance fluctuations and Fano factor) requires some- 
what larger systems. We clearly see in Fig. [2] the relative 
growth of the conductance fluctuations with increasing 
system size and the slow decay of the Fano factor to- 
wards the diffusive 1/3 limit. 



VI. CONCLUSION AND OUTLOOK 

In conclusion, we have analyzed the statistics of trans- 
mission through a sparse chain of tunnel barriers. The 
average spacing of the barriers diverges for a Levy spacing 
distribution p(s) oc l/s 1+a with < a < 1. This causes 
an unusual scaling with system length L (measured from 
the first tunnel barrier) of the moments of transmission 
or conductance, as summarized in Table|TJ A logarithmic 
correction to the power law scaling appears for the first 
moment. Higher moments of conductance all scale with 
the same power law, differing only in the numerical pref- 
actor. As a consequence, sample-to-sample fluctuations 
of the transmission become larger than the average with 
increasing L. 

This theoretical study of a one-dimensional "Levy 
glass" was motivated by a recent optical experiment 
on its three-dimensional analogue [l[. The simplicity 
of a onc-dimcnsional geometry has allowed us to ac- 
count exactly for the correlations between subsequent 
step lengths, which distinguish the random walk through 
the sparse chain of barriers from a Levy walk. We sur- 
mise that step length correlations will play a role in two 
and three dimensional sparse arrays as well, complicat- 
ing a direct application of the theory of Levy walks to 
the experiment. This is one line of investigation for the 
future. 

A second line of investigation is the effect of wave in- 
terference on the transmission of electrons or photons 
through a sparse chain of tunnel barriers. Here we have 
considered the regime of incoherent sequential transmis- 
sion, appropriate for a multi-mode chain with mode- 
mixing or for a single-mode chain with a short coherence 
length. The opposite, phase coherent regime was stud- 
ied in Ref. [1 31 ] - In a single- mode and phase coherent 
chain interference can lead to localization, producing an 
exponential decay of transmission. An investigation of lo- 
calization in this system is of particular interest because 
the sparse chain belongs to the class of disordered sys- 
tems with long-range disorder, to which the usual scaling 
theory of Anderson localization does not apply (22J. 

A third line of investigation concerns the question 
"what is the shot noise of anomalous diffusion" ? Anoma- 
lous diffusion [J] is characterized by a mean square dis- 
placement (x 2 ) oc i 7 with < 7 < 1 (subdiffusion) or 
7 > 1 (superdiffusion). The shot noise for normal dif- 
fusion (7=1) has Fano factor 1/3 pj], 0, and Ref. 
[231 1 concluded that subdiffusion on a fractal also pro- 
duces F = 1/3. Here we found a convergence, albeit a 
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logarithmically slow convergence, to the same 1/3 Fano 
factor for a particular system with superdiffusive dynam- 
ics. We conjecture that F = 1/3 in the entire subballistic 
regime < 7 < 2, with deviations appearing in the bal- 
listic limit 7 — > 2 — but we do not have a general theory 
to support this conjecture. 



Acknowledgments 

We have benefited from discussions with R. Metzler 
and J. Tworzydlo. This research was supported by the 
Dutch Science Foundation NWO/FOM. 



[1] P. Barthelemy , J. Bertolotti, and D. S. Wiersma, Nature 
453, 495 (2008). 

[2] B. B. Mandelbrot, The Fractal Geometry of Nature (Free- 
man, New York, 1983). 

[3] M. Shlesinger, G. Zaslavsky, and U. Frisch, editors, Levy 
Flights and Related Topics in Physics (Springer, Berlin, 
1995). 

[4] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000). 
[5] R. Kutner and Ph. Maass, J. Phys. A 31, 2603 (1998). 
[6] M. Schulz, Phys. Lett. A 298, 105 (2002); M. Schulz and 

P. Reineker, Chem. Phys. 284, 331 (2002). 
[7] P. Levitz, Europhys. Lett. 39, 6593 (1997). 
[8] H. Kohno and H. Yoshida, Solid State Comm. 132, 59 

(2004). 

[9] S. Havlin, A. Bunde, and H. E. Stanley, Phys. Rev. B 
34, 445 (1986). 

[10] K. W. Cheung, K. W. Yu, and P. M. Hui, Phys. Rev. B 
45, 456 (1992). 

[11] E. Barkai, V. Fleurov, and J. Klafter, Phys. Rev. E 61 
1164 (2000). 

[12] The authors of Ref. [TI| calculate a lower bound to 
the mean square displacement, with the result (a; 2 ) > 
£mm(2,3-a) ^ ne m j^j a j position of the particle is ran- 
domly chosen along the chain (so superdiffusion for any 
< a < 2). If the particle starts at a barrier (which 
corresponds to the situation we consider in the present 
work), the result is (x 2 ) > t 2 ~ a (so superdiffusion for 
< a < 1). Earlier papers Qi| gave different results 
for the mean square displacement, but a direct compar- 



ison is problematic because those papers did not notice 

the dependence on the starting position. 
[13] D. Boose and J. M. Luck, J. Phys. A 40, 14045 (2007). 
[14] A. Davis and A. Marshak, in Fractal Frontiers, edited by 

M. M. Novak and T. G. Dewey (World Scientific, 1997). 
[15] H. Larralde, F. Leyvraz, G. Martinez-Mekler, R. Recht- 

man, and S. Ruffo, Phys. Rev. E 58, 4254 (1998). 
[16] S. V. Buldyrev, S. Havlin, A. Ya. Kazakov, M. G. 

E. da Luz, E. P. Raposo, H. E. Stanley, and G. M. 

Viswanathan, Phys. Rev. E 64, 041108 (2001); S. V. 

Buldyrev et al., Physica A 302, 148 (2001). 
[17] C. W. J. Beenakker and M. Biittiker, Phys. Rev. B 46, 

1889 (1992). 

[18] K. E. Nagaev, Phys. Lett. A 169, 103 (1992). 

[19] M. J. M. de Jong and C. W. J. Beenakker, Phys. Rev. B 
51, 16867 (1995). 

[20] We cannot directly take the derivative of Eq. (13. ip . be- 
cause that would lead (for p / 1) to an undefined product 
of 0(L — x) and S(L — x). No such complication arises if 
we take the derivative of Eq. (|3.3p . 

[21] For efficient algorithms to generate random variables 
with a Levy stable distribution, see J. M. Chambers, C. 
L. Mallows, and B. W. Stuck, J. Am. Stat. Ass. 71, 340 
(1976); R. N. Mantegna, Phys. Rev. E 49, 4677 (1994). 

[22] S. Russ, J. W. Kantelhardt, A. Bunde, and S. Havlin, 
Phys. Rev. B 64, 134209 (2001). 

[23] C. W. Groth, J. Tworzydlo, and C. W. J. Beenakker, 
Phys. Rev. Lett. 100, 176804 (2008). 



